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ABSTRACT 


Modifications to the extended Kalman filter and the 
Jazwinski filter are made and compared with the clas- 
sical extended Kalman filter in applications to orbit 
determination using real data. The results obtained 
in this study show that with the kind of data available 
today, the application of filtering theories in this 
field presents many problems. 
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FILTERING THEORY APPLIED 
TO ORBIT DETERMINATION 


INTRODUCTION 

After more than a decade of working with estimation theories and their possible 
applications, no definitive conclusions had been obtained about whether or not 
a filter can be used successfully in the orbit determination field. 

To date most of the studies in this area have dealt with simulated data, and the 
conclusions which have been reached are as follows: 

1. The extended Kalman filter may be used successful^ in the orbit 
determination problem for short lengths of time, i.e., short arc 
lengths. 

2. For applications to longer spans of time, this filter should be modified 
to avoid divergence. This problem occurs because of the lack of 
knowledge of the dynamic system as well as the linearization of the 
original problem. The adaptive filter seems to be the solution to this 
problem. 

This study is an attempt to find a way of applying such estimation theories to 
real world orbit determination problems. 



I 

EXTENDED KALMAN FILTER 

Suppose that the equation of motion is described by the following nonlinear sto- 
chastic differential equation 

dx r = f (x t , r) dr + G (r) dj3 r , r > 0 

Xo~N(i 0 ,P n ) (1_1) 

formally equivalent to: 

E tx 0 l 

E(fx 0 -x 0 Xx 0 -^)'].P o (J 2) 

-^ = f (X T , r) + G (r) u r 

Where: 

x r , the n X 1 state vector 

f (x t , r ) n X 1 continuous function vector 

G (r) , n X r continuous function matrix 

, r X 1 Brownian motion process 
(E jd0 T d/3/ } = Q(r)dr) 

, r XI White gaussian vector process 
(u T ~ N (0, 0(7) ), i.e., E [u T ] = 0, E fu r u/l = Q(r) 

The discrete nonlinear observations can be expressed as 

y K = h(x K )+v K , 

, (1-3) 

v K ~ N(0, o‘ K ). 

(These observations are taken at time instants t K and r = t -t K _ j (r K = t K - t K _ j 
Suppose that we generate a reference deterministic trajectory x t ( satisfying 


{<*,.*><>} 

(u,,r>o) 


1-1 



(1-4) 


d X r 

— I =f(x 7-), T>0. 
dr T 


Define at each point in time 


fix = x 

r r 


— X 


7 


(1-5) 


which satisfies the equation 

d (5x r ) 
dr 


f(x T ,r)-f(x T ,r) + G(r) u r 


(1-6) 


fix Q ~ N (x 0 - x q , P Q ) 

If the deviations (5x) from the reference trajectory are small, a Taylor series 
expansion gives 


d (fix T ) 
dr 


= f 


(x r , r) 5x r + G (r) u t 


(1-7) 


where f x (x r , t) is a function of time only, i.e., (1-7) is a linear system. 


Defining the nominal measurement as 

y K =h(x K ),(x K =x ), (1-8) 

Jv 

then 


5y x =y K -y K = h x ( x k> 5x k +v k>( 5x k =5x r )• 


(1-9) 


Equation (1-7) can be discretized as 


fix K ’ *K - l) ^ X K - 1 + W K 


( 1 - 10 ) 


1-2 



where 4 1 , the state transition matrix, is defined as: 


d » (t, t K _ t ) 
dt 


f x (x T ,T)<i>(t, t K _ J ) 


with 


^(*K- 1 5 *K - l) 7 t_t K-] 


and co K ~ N ( 0, Q K ), where 


i» l K+i 

Q k = $(t K) s)G(s)Q(s)G T (s)^(t K , s )ds 

J *K 


where 

> s ) _ - l) - 1 * s ) 

and 


^(t K _,,s) = $ Hs, t K _,). 

Now the linear filter is directly applicable to the linearized system (1 
1-10) and the equations for 5x and P are as follows: 

a. Prediction: 


^K/K-l - l) - 1 /K - 5 

^K/K-l - ’ *K - 1^ ^ ^ T (*K ’ - l) + 


( 1 - 11 ) 

( 1 - 12 ) 

(1-13) 

(1-14) 

(1-15) 
9 and 

(1-16) 

(1-17) 


1-3 



b. Update: 


6x 


K/K 


= 5x 


K/K - 1 


[ 5y K M X K^ X K/K- l] (1-18) 


P K/K 


1 K t K h x( X K)] P K/K- ! [l — K t R 



(1-19) 


The Kalman gain is 


K. 


l K 


= P 


K/K 


1 h X T (x K ) [' 1 X^ X K^ P K/K - 1 H X T ( x k) + 0 k'. • (1- 


20 ) 


If the reference trajectory is chosen such that 

x(0) = x K ^ 1/K _ j 


( 1 - 21 ) 


we have from (1-5) 


5x 


K- 1/K - 1 


= E { X K - 


1 /K - 1 


"K- 1/K - 1 


-0 


( 1 - 22 ) 


then from (1-16) 


^ X K/K-l 0 


(1-23) 


and 


X K/K- 1 


^| X K/K - E {^ X K/K _ , + x k/K-|| 
X K/K - 1- 


(1-24) 
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In the same way 


P = F 

r K/K - 1 c ' 


and 


|W/K-I X K/K- l)^ X K/K-l X K/K-l) T } 

e {( 5 x k/k-iX 5 x ~k/k-i)} 


e |( x k/k x k/k)( x k/k x k/k) T } 

“ P K/K “ E |(^ x K/K ^ X K/K^} - 


(1-25) 


(1-26) 


Now we may state the evolution of the mean and the covariance matrix of the 
state for the extended Kalman filter: 

a. Prediction: 

x k / K - l obtained from 


dx 

— =f(x r ,r) (1-27) 

dr 

with 

x 0 = X K - 1 /K - 1 (1-28) 

for 


7 T K - l K - 1 


(1-29) 


E K/K-l - > t K- J E K - 1/K - 1 *K - l) + ( 13 °) 


1-5 



b. Update: 


*K/K *K/K-1 + ^t K [^K h(5c K / K _j) 


(1-31) 


The equations (1-19) and (1-20) remain unchanged. 


1-fi 



II 


JAZWINSKI FILTER 

In the orbit determination problem the equations of motion may be stated as 
follows : 


— = f, (R, V, t ) + f 2 (R, V, t) + G(r) u 
dt 1 z 


dR 

dt 


= V 


where: 


fj, includes the accelerations that are well known 
f 2 , possible unknown accelerations and the errors in f j 
On short lengths of time, f 2 can be taken as: 


df 2 

f 2 f 2(r = 0) + dt 


T. 


(r = 0) 


If 


U = f • 


2(r=0) 


df. 


u = 


dt 


(t=0) 


then the equation (2-1) becomes: 


= f, (R, V, r) + u + ur 
dt 1 


if the state noise is not included. 


(2-1) 

(2-2) 


(2-3) 


(2-4) 


(2-5) 


(2-6) 
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The lack of knowledge of u and u makes it necessary to include these in the state 


and estimate them as well as the usual state vector (position R and velocity V). 

Now we may state the complete dynamical system 


— = V 
dr 

(2-7) 

dV 


— — 51 f, (R, V, r) + u + ur 
dr 1 

(2-8) 

du 

— = u 
dr 

(2-9) 

du 

= 0 

dr 

(2-10) 

The initial conditions for r = 0, ( t - t R , ), are 


e{ rT: v t } - e{* t } = 

(2-11) 

E { u } = °K-i 

(2-12) 

e { :, } = 6 k-i 

(2-13) 

E {(x - x)(x - x) T } = P K _ l/K _ l 

(2-14) 

4-#-o) i } = c . k _ 1/k 

(2-15) 

E{(x- S )(n-t,)T^ = C; k i/k t 

(2-16) 

e{(u - a)( u - o) T | = u„ Uk _ i;ki 

(2-17) 
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F {( U ■“)("■“) T l “ U »u K . 
e{(u ■ u)(j • i,) T } =U ; „ K 


1 / K - 1 


1/K - I 


( 2 - 18 ) 


( 2 - 19 ) 


and the measurement model is 


h (x K 5 r K ) + v R , t tj 


( 2 - 20 ) 


where t k is a zero mean, white, gaussian noise sequence with 


e { v k v k T } r k 


( 2 - 21 ) 


The predictions for the estimate of the state (x K/K _ lf u K/K _ u K / K _ and the 
estimate of the state covariance matrix for t = t K , (r = r K = t R — fc K - i ) are: 


K /K - I 


F(X K - 1/K - 1» - 1/K - P “lC/K • 1» t K^ 


U K/K - 1 _ U K - 1 + T K U K/K - 1 


( 2 - 23 ) 


U K/K - 1 U K - i/K - 1 


( 2 - 24 ) 


P K/K - 1 ^ P K - 1/K -1 + <f> C 


1 \p T + ]p C T <f> T 

U K - ] IK - 1 K - 1 /K - 1 


+ \pv 


+ ^U„ 


d/V +<J>C ° \j/ d T + d / d C T • ■ 

UU K-1/K-1 U K ■ 1/K - 1 V 11 K - 1 / K - I 


K ■ 1/K 


l K - 1/K - 


i//d T + ^d u T uu 


K - 1/K - 1 


^ T + ^d U,; * 


( 2 - 25 ) 

UU K - 1/K - 1 


K/K - 1 


= «1>C. 


K - 1/K - 1 


+ ^ u„ 


‘K - 1/K - 1 


r K C- 

K l u K 


+ \p U ■ + \j/ n- • \ 

1/K - 1 UU K - 1/K - 1 d UI, K - I/K - 1/ 


( 2 - 26 ) 
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Cu 


K/K - I ^ ^ u K-1/K-1 + ^ ^ UU K - 1 / K • 1 


+ *d u- 


*K - 1/K - 1 


(2-27) 


where : 


U = U 

UU K/K - 1 UU K - I /K 


+ ^K U uu 


, + t k (U uu 

) 


K - 1/K - 1 


+ U T . 


'K - l/K - 1 


K - 1/K - 


(2-28) 


u • = u • +r v u- • 

UU K/K - 1 UU K - 1/K • 1 K UU K - 1/K - 1 


(2-29) 


u- 


= u- 


UU K/K - 1 UU K - 1/K - 1 


(2-30) 


3x, 


<*> = <!> 


K/K -.1 


3x 


K - 1 


(2-31) 


$ = $ 


d\. 


K/K - 1 


3u 


K -11 


K [ 


\r v i- 


(2-32) 


^ d= ^ d K/K - 1 





4> is the state transition matrix and I, 3 X 3, Identity matrix. 
The update at observation K is: 


X K/K X K/K - 

. + K X 

r" \ 

K < 

X 

i 

( X K/K - 1 T k) 

U K/K = U K/K - 1 

+ K 

s: 

i 

* 

>» 

■ i 

1 1 

< X 

U K/K = U K/K - 


y K M 

X K/K - 1 ’ T k) 


(2-33) 


(2-34) 

(2-35) 

(2-36) 
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P K/K 

P K/K-1 ' 

K x h x Pr/k-i 

c„ 

= c u 

- K. K c „ 

U K/K 

U K/K-J 

^ x U K/K-1 

c, , 

= c- 

-K h C 

U K/K 

U K/K-1 

x x U K/K-1 

u uu 

= u„ „ 

- K K c u 

UU K/K 

UU K/K-! 

l " x U K/K -1 

U-- 

= u- • 

-K, h C- 

K/K 

UU K/K-1 

U X U K/K-1 

u- • 

UU K/K 

= u- • 

UU K/K-1 

h x C u K/K ., 


where : 


(2-37) 

(2-38) 

(2-39) 

(2-40) 

(2-41) 

(2-42) 


9h 


hx hx K UJ X X K/K-1- T r K 

(2-43) 

= = P K/K-i h T x , Y V 

(2-44) 

K u=Ku k =C T u k/k . ] h T x Y-' K 

(2-45) 

K,1=K- k =C t u k/k ., h T x y« K 

(2-46) 

Y K = h x P K/K-1 + 

(2-47) 

The following initial conditions are to be specified 


X 0/0> P 0/0’ U uu 0/0 > U uu Q/0 

(2-48) 

whereas 


A t 

U 0/0’ u o/o 

(2-49) 


are externally set to zero for lack of better Information. 
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The inclusion of Equation (2-42) is the main difference between the filter de- 
scribed here and that of Reference 1. This inclusion avoids the covariance 
matrix of the state from becoming non-positive definite. 
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m 

MODIFICATIONS 


Several modifications have been made to the previous filters in an attempt tp 
solve some problems associated with their performance. 

The main problem related to the EKFis the modeling of state noise. The as- 
sumption of white gaussian noise seems to be poor. But if we accept this 
assumption, we must still assign proper values for the variance. 

On the other hand, the errors of modeling the dynamical system should be 
taken into account to avoid the filter divergence. The optimistic performance 
of the EKF without state noise leads to a growing difference between the actual 
state and the estimated state not reflected in the covariance matrix which 
continues to decrease as more observations are given. The filter divergence 
manifests itself in a secular growth of the predicted residuals which become 
clearly biased. 

In an attempt to avoidthe divergence in the EKF, while maintaining its simplicity, 
a modification has been made in the definition of the Kalman gain. This 
modified filter is going to be called MEKF. 

The definition for the Kalman gain in the MEKF is 


(tj/) K (t K ) (1 + Q K ) 


(3-1) 


where: 


Qk 



+ 



(3-2) 


and 



J_ 

a K 2 r 
0 


if S r j^y K -h(x 
otherwise 


K/K-l >) <0 


(3-3) 
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0 


where: 


" n 2 r-l + n fyK'^K/K-ljJ 


(3-4) 


and K t _ is given in (1-20). 2 is different for each observation type and station, 
and S 0 K = 0. 

This modification tries to avoid the increase in the mean of the residuals 
which are the main manifestation of the divergence. In this way the direction 
is kept but its magnitude is modified to avoid the filter divergence. 

The third filter tested in this study has been a modification of the Jazwinski 
Filter described in Chapter n. 

This filter is going to be called MJF and tries to make the JF adaptive. 

The Jazwinski Filter shows great sensitivity to the values assigned to the 
initial covariances of u and u. It seems desirable to have the ability to set 
these values correctly and automatically. 

The covariance matrix of u can be used to regulate the process, so that the 
residuals are consistent with their statistics. 

Assume that in each step we know U” except for a factor. That is, if U"> is 

UU r J U II 

the covariance matrix given in the process described in Chapter II, the actual 
covariance matrix, U' uu , is going to be: 


U'- • 

UU K-1/K-1 


+ \_ UU 


)Uu„ 

k K-l/K-1 


(3-5) 


and the covariance matrix of the state is 

P k/k-1 = p k./k-i + ^d^uii ^d 1 


We now define 


3-2 



(3-6) 


i [V h (*K/K-i )] 2 

^ + — ■ 

f-1 n ^x P K/K-1 h x T + °k] 


Then 

(h x P K/K-1 h x T+a K) 2 r ==E l r Kp 

where 

r k = [vk ' h ^K/K-l |j’ 

and 




K/K-t "X 


h X +<J K +q K 


h x*d P *d 


T h T 
1L x 


(3-7) 


(3-8) 


(3-9) 


% = 


(^r-D^x P K/K-1 h x T + a K> : 


h x P ^d T h l 


, if positive 


(3-10) 


, otherwise 




n-l 

n 


%-i 



(3-11) 


To avoid the accumulating corrections in successive observations, after 
setting the new value u^ (i ,q K is reset to 


n-l 1 


1- 


n n 




1 +0K 

1 + 8 k 
0 


, if positive 
, otherwise 


(3-12) 
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A further step to completely automate the process would be the initialization 
of If and U — . In the filter tested in this paper this is done by setting 


where K < < 


U„ 


‘o/o 


= Kf, P, 


1 0/0 1 , 


( 3 - 13 ) 




0/0 





( 3 - 14 ) 


1, usually K= 1(T 6 . 
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IV 


COMPUTATIONAL METHODS 


The propagation of the state has been made by numerical integration of the 
Cowell Equations of motion and the state transition matrix by numerical inte- 
gration of the variational equations. 

To obtain accurate prorogation of the state using the least possible computer 
time, two different integration methods as well as two force models have been 
used. 

These two options are: 

1. Numerical integration by a fourth order Runge-Kutta of the equations 
of motion and variational equations with a gravity field including only 
J 2 for equations of motion and variational equations. 

2. Numerical integration by a 12th order summed Adams- Cowell 
predictor-pseudo corrector of the equations of motion and using 
an 8th order summed Adams-Cowell corrector only for varia- 
tional equations, the gravity field being 4 by 4 for equations of mo- 
tion and J, only for variational equations. 

The decision to use one or the other was made by evaluating the number of 
times that the forces had to be evaluated in each case to perform a given 
propagation taking into account the different force model used in each case. 

It was decided to use option 1 only if the number of force evaluations was less 
than or equal to forty for the propagation to the next data point. 

It is possible to have a good estimate of this number by assuming the step size 
used in option (1) satisfies: 


h = QR 2 /m, q = 


0.05 e < 0.2 
0.01 e > 0.2 
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V 


CASES STUDIED 


Three different satellites have been studied in this paper. Two of them have 
very high eccentricity (RAE-B and IMP-J, transfer orbit), and the third is a 
near circular, almost polar orbit (ERTS-1). 

1) The initial conditions for each of these satellites follow. 


a) RAE-B 


Epoch: June 10, 1973 - 

15 hr. 22 min. 

19.9461 sec. 


Initial state: 

a 0 = 198835.0 

e 0 = .9669928 

Km 

x o - “ 

19587.57 Km 


Y 0 = - 

2927.08 Km 

i 0 = 29.11926 

deg. 

z o - - 

8341.75 Km 

£2 0 = 319.3744 

co 0 = 118.5193 

deg. 

x o = - 

3.584823 Km/sec. 

deg. 

Y o = - 

3.735018 Km/sec. 

Mg = 1.321811 

deg. 

z o = - 

2.879267 Km/sec. 


Coordinate system: Mean equator of 1950 


Four different initial covariance matrices have been used. P 03 and P 04 are 
chosen arbitrarily to be diagonal. Pqj and Pq 2 are computed by means of a 
simulation starting with Po ? and Po 4 . The procedure is to start with the epoch 
June 10, 1973, 25 hr 22 min 19.9461 sec and integrate the state backwards 
to June 10, 1973, 14 hr 28 min 47.966 sec. At this time initialize the filter 
with P 0 and P 0 . Propagate these to the original epoch. These propagated 
values are then used for P 0l and P 02 . This procedure helps to eliminate the 
transient response of the filter. 


b) IMP-J 

Epoch: September 26, 1973 - 2 hr. 
Initial state: 


a o 

e 0 

= 124510.0 

= .94719 

Km 

i 0 

= .2867085 

deg. 

n 0 

= 249.7624 

deg. 

^0 

b 120,1121 

deg. 

M 0 

« ,08634 

deg. 


38 min. 15 sec. 


x „ = 

5649.6 

Km 

Y 0 = 

2453.8 

Km 

Z o = 
0 

2434.5 

Km 

2.07466 

Km/ sec . 

10.1934 

Km/sec. 

2.9925 

Km/sec. 


5-1 



Coordinate system: Mean equator of 1950 

The initial covariance matrix (Table 2) was calculated by the Monte Carlo 


method. 





c) ERTS-1 





Epoch: 

August 9, 1972 - 15 hr. 15 min. 0 

sec. 



Initial state: 




a o = 

7283.213 Km 

x o = 

1906.134 

Km 

e 0 

.001863498 

Y 0 = 

-4009.233 

Km 

i-o = 

98.99210 deg. 

Z 0 = 

-5789.823 

Km 

= 

283.5776 deg. 

X 0 = 

.7200720 

Km/sec. 

w 0 = 

148.0186 deg. 

Y 0 = 

5.932595 

Km/ sec. 

M 0 = 

158.4419 deg. 

Z 0 = 

4.338789 

Km/sec. 

Coordinate system: Mean equator of 1950 




Two different covariance matrices (P Q ( , P Q 

) have been used for this case. 


P 0l has been calculated in the same way as Poj and P 0? for RAE-B. The 

state was projected back to August 9, 1973 - 13 hr. 31 min. 0 sec. 

I* 0 has the same diagonal as p but the off- diagonal elements are zero. 

2 05 

2) The observations used for each satellite follow. 

a) RAE-B 

Observation span: 

Start time: June 10, 1973 - 15 hr. 22 min. 20 sec. 

End time : June 11, 1973 - 2 hr. 38 min. 46.254 sec. 

Total number of observations: 233 
Data types: VHF range (p), range rate (p) and 
minitrack direction cosines ( C ,M) 

From a statistical point of view these observations are not good and it is 
difficult to get a D.C, solution from them. Some biases appeared as: 

• CARVON range bias of 480 m. 
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• CARVON range rate bias of 1.32 m/sec. 

• MADGAR range rate bias of .32 m/sec. 

The distribution of these observations is nonuniform. 

The measurement standard deviations are: 

a p = 150 m Oy = 0.0003 radians 

0 - = 20 cm/sec. a = 0.0003 radians 

p M 

b) IMP-J 

Observation span: 

Start time: Sep. 26, 1973 - 3 hr. 2 min. 21.052 sec. 

End time : Sep. 26, 1973 - 11 hr. 35 min. 37.254 sec. 

Total number of observations: 387 
Data type: VHF and Mini track 

These observations are similar in quality to the previous ones. Biases appeared 
smaller and were not taken into account. 

Two sets of measurement standard deviations were assumed: 

1. New: 


0 n 

= 100 m 

o £ = 0.0001 radians 

p 

= 20 cm/sec. 

a = 0.0063 radians 

M 

Standard: 


a 

= 500m 

o E = 0.0003 radians 

p 

= 30 cm/sec. 

a = 0.0003 radians 

M 


c) ERTS-1 

Observation span: 

Start time: Aug. 9, 1973 - 15 hr. 15 min. 0 sec. 
End time : Aug. 9, 1973 - 23 hr. 26 min. 0 sec. 

Total number of observations: 863 
Data type: USB range and angles 
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These observations are the best of the three cases. No biases appeared, but 
the Y angle measurement did not improve the solution. 


Two sets of measurement standard deviations were assumed: 


1. New: 

o X30 = 100 sec. of arc. 

o = 100 sec. of arc. 

y 30 

o X30 = 200 sec. of arc. 

o = 200 sec. of arc. 
y 30 

In all the cases measurements were corrected by 

• Refraction 

• Transponder delay 

• Light time 

The observations were edited by: 

• Elevation angle 

• Amount of the refraction correction 

• k o criterion, i.e., the residual must satisfy 


a p = 10 m 


2. Standard: 


a p = 20 m 


Yk- h (\/ k -i)<K ( h x Ph I +0 k)' y 
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VI 


RESULTS 


The filter errors plotted in Figures 1 through 25 are the result of comparing 
the filter estimate with the DC (iterated weighted least squares) solution. 
These DC solutions have been chosen as "true solutions" for lack of better 
information. The observation spans used for these solutions were longer than 
those used in the filter runs. 

The ERTS-1 DC solution is more reliable than those for RAE-B and IMP-J, 
but the wider scale errors for these two satellites make the results equally 
reliable. In Figures 1 through 20 are plotted differences between the filter 
estimates and the DC solution 

a. Along track: 

(R f -R DC )* V DC I V DC I " 1 

b. Cross track: 

I I R f - R dc I-(Rf - Rdc)- v dc I v dc I"' I 

c. Velocity 

I V F - v DC I 

In Figures 21 through 25 are plotted normalized position and velocity error, 
that is: 


a. Normalized position errors: 

I % -Rdc I - (Pm + P 2 2 + P 33 )" ,a 

b. Normalized velocity errors: 

i v F - v DC i (p 44 +p 55 + p 66 r ,/a 

Figures 1 through 10, and 21 through 25 refer to ERTS-1; Figures 11 through 
15 to RAE-B, and Figures 16 through 20 to IMP-J. 

Figures 1, 2 and 3 compare the errors for the three different filters tested in 
this study, for ERTS-1. The improvement of MEKF appears slight, but both 
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of them show great differences from MJF. From observation 243 to 244 there 
is a time difference of 5 hr. 15 min.; from 516 to 517, 1 hr. 26 min,, and from 
736 to 737, 39 min. These three gaps may be seen clearly in those figures as 
discontinuities. The three filters are sensitive to them but this sensitivity is 
much greater for MJF. This could be expected since the approximation made 
in MJF is good only for short spans of time. For these relatively large gaps the 
covariance matrix of the state grows too much (the factor is ~ t 3 for some 
terms) making the corrections too large. 

Figures 1 and 2 differ in the values assumed for the standard deviations of the 
observations, o p = 10 m, Ox 3Q = ffy 30 = 100 sec. of arc for Figure 1, and 
o p - 20 m,crx 30 =Oy 30 = 200 sec. of arc for Figure 2 (and 2). K is 3 for these 
figures and 100 for Figure 3. The sensitivity due to these is not easily ex- 
plainable unless the characteristics of the observations are taken into account. 

In this case the Y-angle measurement is very poor while range is the most 
reliable measurement. By increasing the standard deviations of all of them by 
the same ratio the filter is opened, this means that more observations are going 
to be included, but this only favors the inclusion of more Y-angle measurements 
causing the relative weight of the range to decrease in some stages of the 
process, thus impairing the solution. 

The irregular influence of this factor may also be seen in Figure 7 (EKF) and 
Figure 8 (MJF). 

The precision reached is about 500 m in position and 50 cm/sec. in velocity 
for EKF and MEKF and about 3 km and 2 m/sec. for MJF, after the first good 
data (range) of the second station are processed. This happens about the 300th 
observation. The precision for EKF and MEKF is intermediate between the 
level of precision used for Range (10 m) and angle (100 sec. of arc ~ 1 km). 

The influence of poor angle measurements seems to be stronger in the MJF. 
u Ji grows too much to adapt to the residuals of this measurement and this is 
reflected particularly in errors in velocity where there is no direct observa- 
bility, coupling in the covariance matrix caused by U" is more difficult to 
control in the velocity states. 

For RAE-B (Figure 11) differences between the filters seems to be smaller. 

The noise used for these observations are o p = 100 m, a- =20 cm/sec., 

°Sl = °m = 0.0003. The initial state is quite close to the actual one. Some kind 
of divergence happens in this case, the position error grows almost steadily 
from 2 km to about 30 km in the 15 hours of process while the velocity error 
grows to 750 cm/sec. at the beginning, and ending at about the same level that 


6-2 



it began, 60 cm/sec. However this kind of divergence is different from the 
result of modeling. In this case it seems to be caused by too few observations. 

For IMP-J (Figures 16 and 17) the noises have been the same as those used 
for RAE-B. In Figure 16, k CT =100 and Figure 17, =10. The results 

shown in Figure 17 are similar to those for RAE-B. The position errors that 
start at 50 km decrease at the beginning as many good observations are proc- 
essed, and grow thereafter to about 250 km at the end of 8 hours of data 
as the observation rate decreases. The velocity errors that started at about 
18 m/sec. grow sharply at the beginning then decrease to a level of 16/sec. 
at the end. 

The results with k CT = 100 (Figure 16) are quite different. After processing 
200 observations and tracking the state better than with k ff = 10 (the position 
error is now 10 km and the velocity error about 1 m/sec.), the performances 
of the different filters change completely. The MJF with a covariance matrix 
which has grown is very sensitive to bad observations and large oscillations 
begin, thus tracking the state very poorly. The EKF and MEKF continue tracking 
the state with approximately the same level of error. In this case the EKF 
shows itself far superior to the other two filters. 

The influence of K 0 seems to be different for each case as the characteristics 
of the observations change. (See Figures 4 (EKF), 5 (MEKF) and 6 (MJF) for 
ERTS-1; Figure 12 (EKF) for RAE-B, and Figures 18 (EKF) and 19 (MJF) for 
IMP-J.) When the assumption of white noise on observations is close to reality 
and enough observations are available as it happens in ERTS-1 data, K a =3 is the 
the best possible value. The inclusion of more observations, or increasing 
the noise on the observations, or increasing K„ deteriorates the filter perfor- 
mances. For RAE-B and IMP-J the assumption of white noise is unrealistic, 
at the same time the observation rate is very low'. In this case the EKF and 
MEKF work much better than MJF and their performances improve from 
= 10 to k„ = 150. For RAE-B the divergence that appeared with k CT = 10 
disappears with k fl = 150, though the error continues being large (10 km in 
position and 80 cm/sec. in velocity). In the IMP-J case the improvement is even 
greater. 

The influence of the initial covariance matrix can be seen in Figures 9 (EKF) 
and 10 (MJF) for ERTS-1 and Figure 13 (EKF) for RAE-B. The EKF and 
MEKF show greater influence to this parameter than MJF does. The adaptive- 
ness of the MJF causes the initial values to be quickly forgotten. This inde- 
pendence is highly desirable since the initial conditions are not well known. 
Figure 13 shows how the covariance matrix calculated as described in Chapter 
V eliminate the transition period. 
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Although the estimation of parameters other than the state is possible and the 
filters tested in this study are being implemented to have this capability, only 
the estimation of biases is available so far. Figure 14 (EKF) and Figure 15 
(MJF) for RAE-B and Figure 20 (EKF) for IMP-J show the difference between 
estimation and no estimation of biases, and it is very simple to deduce that 
there is not enough information in the observations to estimate these other 
parameters. The improvement made by the inclusion of estimation of biases 
is questionable if the data rate is not adequate. 

Figures 21 through 25 show how the covariance matrix reflects the actual errors 
for the different filters in the ERTS-1 case. The errors are normalized with 
their standard deviations as given by the covariance matrix. Assuming that 
filter errors are normally distributed the values plotted in these figures should 
be less than 1 in 66% of the cases. As may be seen in these figures the covari- 
ance matrix is too optimistic reflecting the actual filter errors. The difference 
between actual and expected errors is stronger for EKF, decreases for MEKF 
and is the least in the MJF as may be seen in Figure 21. Figures 22 and 23 
study the influence of k CT . This factor is more important from this point of view 
for EKF than for MJF as may be expected. The assumption of larger observa- 
tion noises makes the covariance matrix reflect a degree of pessimism. 
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(c) 


Figure 3. Comparison of Filter Errors (ERTS-1) 
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Figure 4* 


Comparison of Filter Errors (ERTS-1) 
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Figure 6. Comparison of Filter Errors (ERTS-1) 
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Figure 7. Comparison of Filter Errors (ERTS-1) 
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Figure 8. Comparison of Filter Errors (ERTS-1) 
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Figure 10. Comparison of Filter Errors (ERTS-1) 
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Figure 12. Comparison of Filter Errors (RAE-B) 
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Figure 13. Comparison of Filter Errors (RAE-B) 
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Figure 15. Comparison of Filter Errors (RAE-B) 
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Figure 16. Comparison of Filter Errors (IMP-J) 
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Figure 18. Comparison of Filter Errors (IMP-J) 
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Figure 19. 


Comparison of Filter Errors (IMP-J) 
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Table 1 


Initial Covariance Matrices For RAE-B 











p », 



0.2876E0 

.2284E0 

- .7631E0 

.1207E-3 

.8619E-4 

-.2703 E-3 


.2668 EO 

- .5675E0 

.9615E-4 

.1102E-3 

-.2 144 E-3 



0.2061E1 

-.3191E-3 

- .2088 E-3 

.7275E-3 




.5539E-7 

.4273 E-7 

- .1244E-6 





.5647E-7 

- .9634E-7 






.2966E-6 







.1348 EO 

.1165E0 

- .3458 EO 

.4113E-4 

.2771E-4 

-.8272E-4 


.1508EO 

- .2803 EO 

.3513E-4 

.4302E-4 

-.6686E-4 



.9234 E0 

- .1104E-3 

-.6202E-4 

•2222E-3 




.1631E-7 

.1111E-7 

-.3062E-7 





, 1801E-7 

-.2065E-7 






.6980E-7 




73 

O J 



1. 

0. 

0. 

0. 

0. 

0. 


1. 

0. 

0. 

0. 

0. 



1. 

0. 

0. 

0. 




l.E-4 

0. 

0. 





l.E-4 

0. 

l.E-4 




p », 



0.1 

0. 

0. 

0. 

0. 

0. 


0.1 

0. 

0. 

0. 

0. 



0.1 

0. 

0. 

0. 




l.E-6 

0. 

0. 





l.E-6 

0. 

l.E-6 
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Table 2 


Initial Covariance Matrix For IMP-J 


.889 1E0 -.5432E0 

p », 

.9765 E-l 

.2865 E-2 

-.4460E-3 

.8015 E-6 

.2935E0 

-.6129E-1 

-.1445E-2 

.44 13 E-3 

-.1972 E-3 


.9042 E0 

.2204 E-3 

.3 102 E-3 

.3078E-2 



.28 04 E-2 

•2736E-3 

.1986E-3 

-.6640E-4 
.608 2 E-3 
.2693E-2 


Table 3 


Initial Covariance Matrix For ERTS-1 


.18 17E-2 .2447 E-3 

.4235 E-3 

-.3340E-5 

-.7459E-6 

-.4287E-6 

.2877E-3 

- .2601E-3 

-.4 98 9 E-6 

-.2645 E-6 

-.4263E-6 


.5 127 E-3 

- .7031E-6 

.3132E-7 

.3739E-6 



.6177E-8 

.1408 E-8 

.8776E-9 




.4132E-9 

.4116E-9 
.6443 E-9 
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VII 


CONCLUSIONS AND RECOMMENDATIONS 


Results show that the application of sequential filtering theory to orbit deter- 
mination with the kind of data available today presents many problems. The 
precision reached, although in some cases comparable with that obtained in an 
iterative batch process (DC), is almost always lower. 

The observations have to fulfill some necessary conditions for successful appli- 
cation of the filtering theory to this field, and this condition seems to be more 
strict than those required for a batch process. In this case the problem of esti- 
mation of the state is completely separated from that of estimation of the uncer- 
tainties of that state. In the sequential process these two problems are totally 
coupled so that the real number of parameters to be estimated is bigger in this 
case. The necessary condition that the observations have to accomplish is to 
restore the deterioration suffered by the state and its covariance matrix between 
observations. With an observation rate lower than this, the filter may diverge 
and only with a higher rate than that minimum one, may the estimation of the 
state improve. That minimum rate depends on the orbit type and observation 
type. 

With these conditions accomplished, and with an initial covariance matrix 
perhaps improved b3 r simulating, it seems feasible to use the filtering theory as 
an early orbit determination system in which raw data are used. In this work 
only the USB data used for ERTS-1 were raw. The improvement of the estima- 
tion by using raw data instead of smooth data has not been evaluated. Neverthe- 
less, since in every smoothing process some information is lost, the results 
are going to be improved if this information is retained. 

The potential recovery capability of the MJF has not been tested in this study 
since the possibility of doing so with the data used in this work seemed dim. 

The errors in the estimation of the state make the values obtained for the un- 
modeled accelerations completely unworthy. If the filter does not track the 
state correctly, this means that the unmodeled accelerations are not being 
tracked either, since the two problems are totally correlated. Potential re- 
covery would require a much better tracking system, with almost continuous 
coverage. 

The assumption of white noise for the observation noise is not real in most of 
the cases. This fact has more influence in the sequential processes than in the 
batch process. Another, more exact error model should be found in order that 
successful results in the high precision fast orbit determination problem can be 
realized. 
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